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BPI Observer-based Fault Tolerant Control Design 
for a Class of Fuzzy Bilinear Systems 

Dhikra Saoudi 


Abstract — In this paper, a new methodology of Bilinear 
Proportional Integral (BPI) observer-based Fault Tolerant 
Control (FTC) for a class of nonlinear systems with bilinear 
terms is proposed. Nonlinear systems are represented by the 
Takagi-Sugeno (T-S) Fuzzy Bilinear Model (FBM). The 
objective is to design a stable and robust fault- tolerant 
controller based on a T-S fuzzy BPI observer. The proposed 
BPI observer estimates simultaneously the states of the system 
and the occurring time varying faults with the final goal to 
adapt control law strategies in order that the system remains 
stable even in fault case. Further, the stability conditions of the 
designed BPI observer and the state feedback control are 
analysed with Lyapunov theory to ease the design of FTC. The 
gains of the BPI observer and the state feedback control law are 
obtained by solving linear matrix inequalities (LMI) under 
equality constraints. Moreover, we proposed some sufficient 
conditions for robust BPI observer-based fault tolerant control 
of the FBM with parametric uncertainties. Finally, an academic 
system is presented to illustrate the robustness of the proposed 
controller with respect to uncertainties. 

Index Terms — T-S fuzzy bilinear model, Bilinear 
proportional integral observer, Observer-based control, Fault 
tolerant control, Lyapunov function, Parametric uncertainty, 
LMI. 


controller designs based on fuzzy observers for T-S fuzzy 
systems. Nonetheless, in [17], [18], the authors have 
proposed sufficient design conditions for robust stabilization 
of T-S fuzzy models with parametric uncertainties based on 
state estimation. However, to the best of our knowledge, no 
previous study has investigated the problem of an BPI 
observer-based FTC for T-S fuzzy bilinear models and a 
robust stabilization for uncertain T-S fuzzy models. 

In this paper, we consider the fault-tolerant control 
problem for the fuzzy bilinear systems with parametric 
uncertainties. Fuzzy BPI observer is designed to estimate 
simultaneously the states and the occurring time varying 
faults. Based on the fault estimation, the fault-tolerant 
controller is designed to ensure that the closed-loop system 
with faults is asymptotically stable. 

This paper is organized as follows. In section II, the 
considered structure of the T-S fuzzy bilinear system is 
presented. In Section III, the problem statement is described. 
In Section IV, the problem of observer-based control for the 
FBM is developed. In section V, the robust control of the 
FBM with parametric uncertainties is proposed. An example 
is provided to show the effectiveness of the proposed design 
in the subject of section VI. 


I. Introduction 

The Bilinear systems are a special class of nonlinear 
systems, where the dynamics are jointly linear in the state and 
input variables, that attract the interest of researchers [1], [2]. 
Moreover, there are many engineering applications as well as 
models in biology, ecology, nuclear engineering are naturally 
described by bilinear systems. Furthermore, many nonlinear 
systems can be adequately approximated bilinear systems. 

As a special extension, fuzzy bilinear systems based on the 
T-S fuzzy model with bilinear rule consequence have 
attracted the interest of researchers [3]-[5]. It is proved in 
these papers that often nonlinear behaviors can be 
approximated by T-S bilinear multimodel description. Some 
topics of control are extended to T-S FBM, such as stability 
and stabilization in [6], [7], observers and state estimation in 
[8]-[10], faults detection and fault tolerant control [9]-[ll]. 

Fault Tolerant Control have become challenging problems 
in the area of modern control theory. FTC allows having a 
control loop that fulfils its objectives when faults appear. 
Thus great interest has been accorded to the fault tolerant 
controllers for nonlinear systems in the literature [12] -[14]. 
Moreover, it is noted that all of the aforementioned assume 
that the system states are measured, which is not true in many 
control systems and real applications. For this reason, 
observer-based fuzzy controllers were considered in many 
researches. In [15], [16], the authors have studied the 
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II. T-S Fuzzy Bilinear Systems modeling 

The T-S fuzzy bilinear models based on the T-S fuzzy 
model with bilinear rule consequence are defined by 
extending the T-S fuzzy ordinary model. Similar to [19], the 
fuzzy bilinear model can be represented by the following 
fuzzy if-then rules: 

R l : if £(f) is F n and- --and £(f) is F ig 

then f = A x(t) + B.u(t) + N t x(t) u(t) ^ 

[y(0 = Cx(t) 

where R denotes the i th fuzzy rule V/ = {l,« • •, r} , r is the 
number of if-then rules, ( t ) are the premise variables which 
can be measurable or not measurable, and are 

fuzzy set. x(t) e SR n is the state vector, u (t ) e 91 is the 
control input, and y(t) <e91 p is the system output. The 
matrices A i , B j TV , C are known matrices. Then, 
the overall FBM can be described as follows: 


x(t) = h t . (^(0) ( A i x (f) + BjU(t) + N t x(t) u(t)) 

i=\ 

y(t) = Cx(t) 


( 2 ) 


where h.(.) verify the following properties: 


$> (£(*)) = 1 


;= i Vi e{l, 2, •••,?•} 

0 < < 1 


( 3 ) 
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In the remainder of the paper, the following lemmas are used. 

Lemma 1: For any matrices A and B with appropriate 
dimensions, the following property holds for any positive 
scalar X •’ 


p is a given scalar, and D e %l lxn , i = 1, • • • , r are vectors 
to be determined. 

In the following section, we will be developing the BPI 
observer-based control problem for the T-S FBM. 


A t B + B T A < AA t A + A l B T B 


( 4 ) IV. Stabilizing fault tolerant control for FBM 


Lemma 2: Given a scalar p> 0 and a symmetric positive 
definite matrix P , the following inequality holds: 

2 x T y<p 1 x T Px + py T P l y, x, y e 9T (5) 

Lemma 3: Let us consider P a positive definite matrix and 
Q a full column rank matrix. It follows that the matrix 
QPQ t is a positive definite matrix. 


III. Problem statement 


Let us consider the T-S fuzzy bilinear models containing 
the actuator faults which can be rewritten in the following 
form: 


x(t) = (#(*)) 

i = 1 


r A t x(f) + B t u(t) ^ 
K +N t x(t)u(t) + G i f(t)j 


y(t) = Cx(t ) 


( 6 ) 


where f(t) is an actuator fault. 


In order to estimate both state and faults of FBM, we 
propose a BPI observer and make an analysis of the error 
system, from which we provides a design method of a fuzzy 
observer for the fuzzy bilinear system (2). The considered 
BPI observer is described by the following equation: 

H,z(t) + L,y( t) + J i u(tV 


z(0 = &(£«) 


+M i u(t)y(t) + <p i f(t) 


x(t) ■■ 
y(t) 


■ Z(t)-Ey(t ) 

■ Cx(t) 


( 8 ) 


f(t) = yy, (4V))'i (KO- y(0) 


where z (t) e 9T , x(t) e 9T , y(t) e 5R P and f(t) e 5R 7 are 
the observer state vector, the estimated state vector, the 
estimated output vector and the estimated faults vector 
respectively. H i , M i , F , J t , cp t , x ¥ i and E 

are unknown matrices of the BPI observer of appropriate 
dimensions to be determined. 


Assumption 1: The fault f(t) satisfies < a x and the 
first time derivative of f(t)with respect to time is norm 
bounded i.e. \\f(t)\\ < a 2 and 0<cq<oo,0<cr 2 <oo. 


Note that the premise variables do not depend on the state 
variables estimated by a BPI observer. 

Let us consider the state and fault estimation errors defined 
respectively by: 


The main objective of FTC design is to find a control law 
u(t ) such that the closed-loop system remains stable despite 
the presence of actuator faults. Then, it is necessary to 
estimate the states and faults. In this paper, a BPI observer is 
utilized for estimate both the states and the actuator faults 
with the final goal to adapt control law strategies such that the 
system remains stable even in actuator fault case. Following 
the design concept in [1], the fuzzy control law for FBM is 
formulated as follows: 


M(f) = XV£W) 

i=\ 


pD,x(t) 

yjl + x T it)Df D t x(t) 


= y^h l (%(t))psin0 i (t) 

J=1 


(7) 


= XA pD,x{t) cos Oft) 

i = 1 


where: 


sin 6 j (t) 


cos 6 t ( t ) 


DjXjt) 

7l + x r (r)D, r D ; x(0 
1 

^l + x r (OA r £>x(0 


6 j e 


71 71 

~2 ’ ~2 


e x (t) = x{t)- x(t) (9) 

e f(t) = f(t)-f(t) (10) 

From (8) and (9), the state estimation error given as follows: 
e x { t ) = Tx {t)- z{t) (11) 

where p = I n + EC • 


By taking into account (2) and (8), the dynamics of the state 
estimation error and the closed-loop system with the state 
feedback control law defined in (7) becomes after some 
calculations: 


«.(0=2>(£(o) 

i = 1 


/ H i e x (t) + (TA i -LC-HT)x(t) 

+ (TN i -M i C)x{t)u{t) + (TB i -J i )uit) 
+ (TG l -<p i )f(t) + <p i e f (t) 


( 12 ) 


r r 


i = 1 7=1 


dT.xO) - pBPj cos 6 j e x ( t 
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where: 

dx. = A + pN\ sin 6j + pB j D j cos 6 j 


(13) 


(14) 


Then, the state x(t) of the system, the state estimation error 
e x (t) and the fault estimation error e (f) are bounded. The 

gains of the observer and the state feedback control law are 
given by: 


Hence, if the following conditions hold true 
Vi e {l, 


E = P~ l S , ./ = P~ 1 R i , M, = P~ l U l , 

V = (L + EC) A, - P 2 'W,C , p, = P~% , 

L l = p;'w : - H t E , x ¥ j = p;% , D [ = VtQT 1 


(25) 


TA i -LC- HT = 0 
77V, —M;C = 0 


TB i -J i = 0 


(15) 

(16) 


Then, the equation of the observing error becomes 

4 ( t ) = Y J h l {%(t)){H i e x (t) + , Z i f(t) + <p i e f (tj) 

1=1 

with s . = TG i - cp t 

From (15) and using T = I + EC , we get 
H i = TA l - K C 
with K, — H.E-\- L. 


Proof: In the framework of a FTC by state-feedback, a new 
procedure is proposed such that the closed-loop system and 
the convergence of the state and fault estimation errors 
becomes asymptotically robust stable by using a Lyapunov 

(17) function-based design approach. Let us consider a Lyapunov 
function V(t ) depending on x (t ) , e x (t) ande f (t) defined 

by: 

(18) V(t) = x T (f)P x x(f) + e T x ( t)P 2 e x ( t ) + e T f ( t)P 3 e f ( t ) 

(26) 

where p i , P 2 and P 3 are symmetric and positive definite 

matrices with appropriate dimensions. Stability condition for 
the estimation error yields that the time derivative of (26) 
should be negative definite. The derivative of V(t) with 

(19) respect to time yields : 


( 20 ) 


The goal is to design a feedback control law such that the 
system remains stable even if a fault occurs. A result is 
summarized in the following theorem. 


i=l 7=1 


X (f)FJx(O + X ( t)P x x(t ) 
+e T x (t) P 2 e x (t)+e T x (t) P 2 e x (; t ) 
+e T f 0 t)P 3 e f (; t ) + e T f ( t)P 3 e f (t) 
(27) 


Theorem 1: If there exist symmetric positive definite that becomes by using equations (18) and (13) 
matrices Q,T l ,T 2 ,T 3 ,P l and P 2 , matrices w t , V t 
Y. , Z. , R i , U i , S and positive scalars s , J3 , p> 0 
such that LMIs (21) subject to linear equality constraints 


22), (23) 

are 

satisfied V/ 

= 1 • • 

• r , 

ii 

>* 

• • r • 

LMl 

* 

* 

* 

* 

* 

* 

* 

N t Q 

-/ 

* 

* 

* 

* 

* 

* 

By j 

0 

-I 

* 

* 

* 

* 

* 

Q 

0 

0 


* 

* 

* 

* 

p^yf 

0 

0 

0 

-20Q 

* 

* 

* 

0 

0 

0 

0 

0 

-2/31 

* 

* 

0 

0 

0 

0 

pi 

0 

LMl 

* 

0 

0 

0 

0 

0 

PI 

LM3 

£ l T 2 


( 21 ) 


<0 


/=1 7=1 


where: 

n,=o^+^ 


n=H'p + p 2 H 


x T n..x + e T x Cl.e x + 2x T P.G.f 

IJ XIX l l j 


-2 px P^B.D. cos 6 e y 
+ P P 2 E-if + 2 e x P 2 <Pi e f 
+2e T f P 3 f + 2e T f P 3 'V,Ce x 


(28) 

(29) 

(30) 


r,=(p 2 +sc)b, 

U i C={P 2 +SC)N i 


(22) By using the Assumption (1) and applying Lemma (2) for 
three terms of the above inequality (28), it follows that: 

(23) 


LM\ = A i Q + QA' i +p 1 I 

LM 2 = ( P 2 + SC ) A +Aj(P 2 + SC f - W t C - (Wf) T + e~ 1 T i 
LM3 = Y, T +Z ; C 


(24) 
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P 2 E,f < e-'eXe x + ef (sf P 2 T; , P 2 E i ) f 
<s-'e T J,e x+ ,j 

2e T f Pj < e~V f r 2 e f +sf T (P£?P 3 )f 
<£~V f r 2 e f +rj 2 

2 X T P X GJ < s - l x T r 3 x + sf T (Gf Pf-'Pfr ) / 
< £~ 1 x T T 3 x + rj 3 

with: 

ri 1 =sa 2 i ^{E T i P 2 T- l P 2 E i ) 

V ~ Sa 2 Ailax (^3^2 ^3 ) 

rf =£afA_{GjP^Pfi i ) 


r = ™o nA '"‘" 


V * =1 7=1 


(31) 


which can be also bounded by: 


t < min 


M- a «) 


If the following inequality holds: 


zi^w^(o)A # <o 

/=1 7=1 

(32) Then, we can obtain that 

V{t)<-T\\X a {t)f +s 


( 38 ) 


(39) 


(40) 


(41) 


Let us define the scalar 8 the maximum value of the sum of 
all constant coefficients obtained in the equalities (32): 

8 - max (77 1 +rj 2 + 77 s ) (33) 

Then, the time derivative of the Lyapunov function (28) is 
bounded as follows: 


i= 1 7=1 


x T Yl ij x + e r x O, i e x 


—2px r P l B i D j cos 6 j e x 
+2e x P 2 < P i e f + iCe x 
+£~ 1 e T x rie x +£~ l e T f r 2 e f 
+£ ^X T P^X + 8 


It follows that v(t) < 0 if r||X a (0f >8 ,\/t> 0 • Then, we 
can deduce that the state x(t ) , the state estimation error 
e x (t) and the fault estimation error e f ( t ) converge to a small 

set according to Lyapunov stability theory and lie in it. This 
set is smaller as the constant 8 converges to zero. 

Let us introduce the following notations for the sake of 
simplicity: 

Y,=tlh l (m)Y, < 42) 

i=l 

Y^-iii^rnpAmY,, (43) 

i= 1 7=1 


The above inequality (34) can be reformulated as follows: 

(35) = 


(34) where Y t and T. are given matrices. By using these notations, 
the inequality (40) becomes: 


Y(t) < X T a (r)Z2> (ZWfo («0)A V x a (t) + S 

i= 1 7=1 


where the augmented vector X and the matrix A., are 


respectively given by: 

x a (t) = 


x(t) 

e x (t) 

e f (t) 


A,, = 


n, r 3 cos 6. 


-ij ■ - - 3 
* 


0 


n,+s~ l r, PM+i'i’.cf p 3 

* ^r'r. 


©' © 2 

* © 4 


<0 


with: 


©k - n y + T 


(36) ©i =[-/3/)fi,D j cos^ 0 


© 

(47) 


+( 4 P,c) r p 3 

* 


Let us consider a symmetric matrix % defined as: 


(44) 


(45) 

(46) 


Let us define the positive scalar r 


(37) 


X = 




0 

Xi 


>Zi = 


p 1 

r i 

0 


0 

I 


(48) 


By using Lemma (3), and post and pre-multiplying the 
inequality (44) by % , we can obtain that: 
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P < ' 0 e^' 0 (49) 

* \ 

The term xft%X\ can re pl ace d by considering the 
following inequality which holds for any scalar (3 such that 


(50) 

(51) 


bilinear model with actuator faults: 




y(t) = Cx(t) 


(A + AA ) x(t) + {B i + AB, ) u(tj 
+ (N t +AN i )x(t)u(t) + G i f(t ) 


(56) 

The parameter uncertainties considered here are 
norm-bounded and presented by the form 

[M A Bi A A] = ti,F,(t)[E li E 2i E 3i ] (57) 


Then from the above inequality (51) and with the Schur 
Complement, it follows that the inequality (49) holds if: 

V © 1 ^ -1 o 

* ~2j3 Xl pi 
L & 

Now, the term of (52) can be rewritten as: 


<0 


(52) 


where %. , E u , E 2i , E 3i are known real constant matrices of 
appropriate dimensions, and F t (t) is an unknown matrix 
function satisfying F T < I , in which I is the 
identity matrix of appropriate dimension. 

The feedback controller and the state of BPI observer for T-S 
FBM with parametric uncertainties (56) is the same as that 
for (6). 


= P~ l Af + A P~ ] + pNf- 1 sin 0. + pP~ l Nf sin 0. 

+ pBjDjPp cos 0 j + pPp (B'D J cosd j +£-'p-'r-.Pp 

(53) 

Applying lemma (1), the previous equation can be rewritten 
as: 

p-'vW Z + A p i~' + P~ l Nf N'P ; 1 + p 2 i (54) 

+ ^{3,0^ bee;' +£ -' p;'r i pp 


In order to describe the dynamic of BPI observer (8), the state 
and fault estimation errors is defined by (9), (10). Then, from 
(56) and (8), we have: 


' H e x (t) + (TA l - pC - H,T)x(t ) 


^(0 = Z A .-(^w) 

1=1 


+ {TN i -M i C)x(t)u(t) + (TB i -J i )u(t) 
+T A A { x{t) + T A B t u(t) + T A N t x(t)u(t) 


+ (TG l -(p,)f(t) + cp i e f (/) 


(58) 


From the inequality (54) and the matrices (46), (47), using the 
notations (42), (43), replacing /f by (19) and T = I n + EC 

and by introducing the following variable changes: 


r r 


i = 1 M 


(^y + A®,y ) X(t) + G t f (t) 

-p ( B. + A B, ) D ] cos 0.e x (0 y 


Q = p;\ Dj = VjQ , W t = P 2 K. , Y, = P 2 ( Pl 
p , = p P , U, = PM, , Z, = PpV, , S = P 2 E 


(55) 


where: 


(59) 


we can obtain the inequalities given in theorem (1) by 
applying the Schur's complement [20] to (52) under equality 
constraints (22), (23). This completes the proof of the 
theorem. 


dv = A. + pN j sin 6j + pB i D j cos 6 } (60) 

AO y = A A. + pAN t sin 6 } + pAB i D j cos 6 } (61) 


This section is devoted to the design of a BPI observer-based 
control for T-S FBM. The system parameters used are known 
but in reality the system parameters may either be uncertain 
or time-dependent. Moreover, the problem of stabilization 
remains a key problem in the study of uncertain T-S fuzzy 
control systems and their controller designs. The study of BPI 
observer-based control for T-S FBM with parametric 
uncertainties will be subject the next section. 


If the conditions (15), (16), (17) are satisfied 

\fi g {l, 2,* • *, r} , we get: 




r H,e x (t) + TAA,x(t) + TAB,u(t)"' 
+T AN t x(t)u(t) 

+ {jGi - (p^ f (0 + (p t e f (t) 


(62) 


V. Robust control of the fuzzy bilinear model with By using (7), we obtain: 

UNCERTAINTIES 

In this section, we shall consider a uncertain nonlinear 
system which can be described by the following T-S fuzzy 
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r r 


i=\ 7=1 


^k(0 + Sf/(0 ' 

s +’E? i x(t) + <p i e f (0 y 


where : 

E)=H i -pTAB i D j cos (9. 

$=7G-n 

Z 3 t = TAA t + pTAB i D j cos £ . + pT AN { sin 0. 


(63) 



M'XK 

* 

* 

* 

* 

* 




E 2 ,Dj 

- M I 

* 

* 

* 

* 




E 3i 

0 

-pi 

* 

* 

* 


LM6 - 

ph ] 

0 

0 

-pi 

* 

* 




E 2i D 2 

0 

0 

0 

~p-'i 

* 

(64) 



_p+2p 2 ) 

n] 0 

0 

0 

0 

-p 7 




~ —2fiX\ * 






(65) 

LM1 = 

LM3 a 








LM2 * 






(66) 

A = 

LM 3 e~'T 2 







Hereafter, we give the following theorem which presents 
sufficient conditions such that the closed-loop fuzzy system 
(59) is globally asymptotically stable in the presence of 
actuator faults and which gives the gains of the observer and 
the state feedback control law. 


Then, the state x(t) of the system, the state estimation error 
e x (t) and the fault estimation error ^ / (^)are bounded. The 

gains of the observer and the state feedback control law are 
given by: 


Theorem 2: If there exist symmetric positive definite 
matrices Q,T l ,T 2 ,T 3 ,P l and P 2 , matrices w t , V t 

Y. , Z. , R t , U t , S and positive scalars z , , p , p> 0 

such that LMIs (67) subject to linear equality constraints 
(68), (69) are satisfied Vz = 1- • • r , V/ « 1- • • r . 


LM\ 


** 


* * * * 


N.q -i * * 

By j o -/ * 

Q 0 0 -T 3 1 

E U Q 0 0 0 

E 7 y j 000 
e 3i q 0 0 0 

LM4 0 0 0 

E 2i Vj 0 0 0 

LM5 0 0 0 

p(B i V j ) 0 0 0 

0 0 0 0 




* * * 

* * * 

* * * 


—Jjj * * * * 

0 -pi * * * 

0 0 -jul * * 

0 0 0 -pi * 

0 0 0 0 -ft-'l 

0 0 0 0 0 -/f'l * * 

00000 0 LM6 * 

0 0 0 0 0 0 0 IM7 



R l ={P 2 +SC)B i 

U l C={P 2 +SC)N l 

where: 


(67) 

( 68 ) 

(69) 


LM\ = A i Q + QA t T + p 2 l 

LM2 = (P 2 + SC)A i + Af (P 2 + SCf -Wf - {W,C) T + ^‘r, 

LM2 = Yj + Z,C 

LMA = ph](l + EC) T 

LM5 = \+2p)h] (/ + ECf 


E = P~ l S , J, = p-'R , M, = P-;'U l , 

H, = + E C) A, - P-'W,C , <p, = P~% , (71) 

L , = P-I% - H t E , T, = P~% , D = V,Q-' 

Proof: In order to prove the stability of the closed-loop 
system and the convergence of the state and fault estimation 
errors, sufficient conditions are derived using Lyapunov 
function (26). Then, the derivative of y(^)with respect to 
time yields: 


i = 1 7=1 


x T (n .. + An,.. ) x + e T x n v e x 
+2x T P l G l f + 2e T x P 2 <p l e f 
-2 px T P 1 ( B i + A B i ) Dj cos 0,(' x 
+2e T x P 2 E 2 f + 2e T x P 2 ^x 
+2e T f Pj + 2e T f P 3 '¥ i Ce x 

(72) 


where: 


I I;/ - 7J ( I>„ 


(73) 


An^AO.+A®^ 

(74) 

Qj = ( Hf — p(TEB i D j f cos ) P 2 
+ P 2 (//, - pTAB i D j cos 0j ) 


(75) 


By using (31), the time derivative of v(t) (72) is bounded as 
follows: 
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;=i j = l 


x T (fU + Ally ) x + e r x Q. ij e x 
-2 px r P l (B i + AB i ) Dj cos 6 j e x 
+2e T x P 2 (p i e f +2e T x P 2 E 3 i x 
+ 2e)P^ i Ce x + 6~ l e x T l e x 

+£* 6 j-P 2 6 j + G X + 


(76) 

where S is defined by (33). 

The above inequality (76) can be reformulated as follows: 

V(t) < X T a ( t)YZ h i (^(0)A „x a (t) + s ( ?? ) 

i = 1 7=1 

where the augmented vector X a is defined by (36) and the 
matrix A. is given by: 


A ,= 


A‘. A 2 0 

* V, Ay 

* * e l r 0 


with 
A 1 . = n 


, Any 


£ -T, 


A* = -pP, (B, + A B, ) D } cos 6j + (S, 3 f P 2 
(80) 

A 3 9 =P 2 ip i +{'¥ i C) T P 3 
From (91), we have: 


A y = Ay + Ay 


where a . is defined by (37) and : 


\ij A 3 ij 


Ay - ^ 

* 

(83) 
with: 


Aw-AHy 


2 ij 
* 


A 2 ij - pP 2 T A BiDj cos 6- + p{T A Bfl . ) P 2 cos Qj 
Kj = -pPAB.Dj cos dj + (Sf ) r P 2 


Ajy. = A A i P l 1 + pAN j P l 1 sin Oj + pAB i D j P l 1 cos Oj 
+P t 'AAj + pP~ l ANj sin Oj + pP t ' (Afl,D ; ) T cos <9, 
A 2 y = pTAR Dp cos 0, + pP 2 -‘ (PAfi,D. )' cos 0 ] 
A 3 y = -pABPjP- 1 cos Oj + P 2 -‘ (PAA ) r 
+pP 2 1 (r ABjDj ) r cos <9 ; + pP 2 ' (TAN i )' sin 0, 

The equation A tj can be reformulated as follows: 

Ay=A‘ + Aj + Aj 

where: 


( 85 ) 


( 86 ) 



0 

A 11 

0 


0 

A 21 

0 


IJ 


y 

A 1 . = 

V 

* 

0 

0 

>Al = 

* 

0 

0 


* 

* 

0 


* 

* 

0 




~A 

uj 0 

o" 




(87) 


(78) 


(79) 


Ay = 


a 2iJ 0 


% = P- 1 (TAA j ) T + pP 2 1 (PA/?D ; ) r cos Oj 
+ pP; , (TAN i f sin 6. 

Ay 1 = -pAB i D j P l ' cos Oj 


( 88 ) 


Using the uncertainties structure and by using the separation 
lemma 1 [21], we obtain 


(81) 


(82) 


Ay = 


V 

* 

* 


0 0 


2 
7 

* 0 


T, 0 


(89) 


(84) 


where: 

=(\ + 2p 1 )pTh l F i Flfi r i T T + M P l - l D]E T 2i E 2i D ] P l - 1 
+p~ 1 Pj 1 Ej i E u Pj 1 + p'P;'D] E T 2i E 2i DjP~ x 
+p ' P~' EjE^Pj' + p 2 p~ 1 Th i F i F l T ti r i T T 
f 2 = p'p-'ElE u p-'+p-'P 2 'ElE 2 ,P 2 ' 

+P'P-'D] E(.E r p.PA + fi-'p'hftft?, 

+ (l + 2 P 2 ) ph i F i F i T h T i + pP 2 x D] E T 2i E 2i DjP 2 ~' 

From (82) and (89), we have: 


(90) 


0 1 © 2 

* © 4 


0 1 0 2 

^CC W cc 

* © 4 


<0 


(91) 


Multiplying the previous equations (84) on the left and right 
by p~ l and p~ l , we can rewrite (84) as: 


where ©| , © 2 ? ©* are defined by (45)-(47) respectively and 
0 1 © 2 © 4 are given by: 

V ’ U ’ V ° 
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©H° o] 


(92) 


Using (49), (51) and with the Schur Complement, it 
follows that the inequality (91) holds if: 


p-W p^+C) 1 


0 

0 

* 


0 

0 

* 

* 

~2PZi 

pi 

* 

* 

* 

0 4 


We calculate each terms of inequality (93) and by 
considering the variable changes (55), we can obtain the 
inequalities given in theorem (2) by applying the Schur's 
complement [20] to (93) under equality constraints (68), (69). 
This completes the proof of the theorem. 


VI. Numerical example 


To illustrate the effectiveness of the proposed method, an 
example is considered. The dynamical FBM is given by the 
following data: 


"-2 

1 

0 " 


" 5 " 


"-2 

0 

0" 


v 

1 

-3 

0 

• a = 

1.2 


0 

0 

0 

, G l = 

0 

2 

1 

-8 


2.5 


0 

0 

1 


i 



"-2 

0 

1 " 


" 4 " 


"-2 

0 

0~ 


v 

4 = 

0.5 

-2 

0 

,b 2 = 

2 

’ N 2 = 

0 

0 

0 

, G l = 

0 


2 

1 

-1 


0 


0 

0 

1 


i 


Uncertainties are also defined by the following matrices: 



"0.1 

0 

0 " 


“0 

0.5 

0 " 


" 1 " 


"0 

0 

0“ 

»i= 

0 

0 

0.1 


0 

0 

0.3 

p - 

’ ^21 

0.5 

p - 

’ jl 

0 

0 

0 


0 

0.1 

0 


0 

0 

0 


0 


0 

0 

0 


"0.1 

0 

0" 


“0 

0.5 

0 " 


~ 1 " 


“0 

0 

0 " 


0 

0 

0.1 

’ E\2 - 

0 

0 

0.6 

f - 

’ ^22 

0.5 

P - 

, ■C'32 

0 

0 

0 


0 

0.1 

0 


0 

0 

0 


0 


0 

0 

0 


The weighting functions are defined by: 


\ (x(t)) = exp 




fy(x(o) 


The additive actuator fault signal is defined as follows : 


0 




0.5 sin (0.5;r) 
0.3 
0 


t<ls 

ls<t<ls 

ls<t<9s 

t>9s 


Thus, the resolution of the conditions of theorem 2 leads to 
the following controller and observer gain matrices 
respectively with p = 0.05 


A = 

[3.343 -0.098 

1.254] , D 2 

= [-0.041 


"9.022 -4.533" 


"9.022 0" 

E = 

5.976 -8.965 

,m, = m 2 = 

5.976 0 


2.046 -3.569 


3.046 0 



' 0.688 

-1.056" 


'-504.851 

ii 

o 

* 

-0.399 

0.584 

. h= 

656.529 


-1.567 

2.344 


165.730 


-1.008 -0.254] 


752.266 

-987.282 

-249.618 



-18.166 

-9.956 

t-H 

T 

j 1 =[-4i.m - 

-27.882 - 

-13.229f 


'-51.288 

-10.218 

-77.669" 


' 23.644 

46.592 

-0.148" 

H t = 

29.304 

31.204 

-32.793 

,h 2 = 

-53.635 

-17.602 

-25.086 


112.873 

80.082 

15.084 


-14.954 

-1.058 

-9.543 


(p x - (p 2 =\ 4.511 2.988 1.523f 
='P 2 =[-1.523 -2.988] 


The simulation results are shown in figures (l)-(4). The 
time-varying actuator fault and its estimate are depicted on 
Figure 1. It appears clearly the good estimation of this 
actuator fault. Furthermore, the BPI observer provides the 
fault estimation which errors are illustrated in the figures 2. 
One can shows the fault estimation error is zero-mean except 
at time Is, 7s and 9s that suits abrupt changes of the fault. 



Fig 1 . Evolution of the fault and its estimate 
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Fig 3. The output y x (t) of the system: nominal output, 
output without FTC and output with FTC 
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Fig 4. The output y 2 {t) of the system: nominal output, 
output without FTC and output with FTC 

The figures (3-4) compares the outputs of the nominal model, 
the outputs of the faulty system without FTC and the outputs 
with FTC. One can see that the trajectory of the system 
follows the trajectory of the reference model, even in fault 
case. It is clear that the proposed strategy is robust with 
respect to the actuator additive fault f(t) and with parametric 
uncertainties. Moreover, the BPI observer-based FTC 
scheme provides good results in the presence of 
norm-bounded parametric uncertainties, as shown in the 
figures. 


VII. Conclusion 

This paper proposes a state feedback fault tolerant control 
law methodology for T-S fuzzy bilinear systems. A T-S fuzzy 
BPI observer is designed for the proposed FTC strategy in 
order to simultaneously estimate time varying faults and state 
variables. This method aims to adapt the control law in order 
to compensate the effect of the time -varying faults by 
maintaining the stability of the system. The stability 
conditions of the designed observer and the state -feedback 
control has been provided and solved through a set of Linear 
Matrix Inequalities under equality constraints. Moreover, 
sufficient conditions for controller design based on state 
estimation for robust stabilization of T-S FBM with 
parametric uncertainties has been proposed. Finally, the 
performance of the proposed BPI observer based FTC 
approach has been given to an example in order to illustrate 
the validity of this method. 
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